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ABSTRACT 

The Smithsonian Hectospec Lensing Survey (SHELS) is a magnitude-limited spec- 
troscopically complete survey for R < 21.0 covering 4 0°. SHELS provides a large 
sample (15,513) of flux-calibrated spectra. The wavelength range covered by the spec- 
tra allows empirical determination of k-corrections for the g and r bands from z = 
to ~ 0.68 and 0.33, respectively, based on large samples of spectra. We approximate 
the k-corrections using only two parameters in a standard way: D„4000 and redshift, 
z. We use Dn^OOO rather than the standard observed galaxy color because L'n4000 is a 
redshift-independent tracer of the stellar population of the galaxy. Our approximations 
for the k-corrections using D„4000 are as good as those based on observed galaxy color 
(9 ~ ''') of the scatter is ^ 0.08 mag). The approximations for the k-corrections are 
available in an online calculator. Our results agree with previously determined ana- 
lytical approximations from single stellar population (SSP) models fitted to multiband 
optical and near- infrared photometry for galaxies with a known redshift. Galaxies with 
the smallest Z)„4000 — the galaxies with the youngest stellar populations — are always 
attenuated and/or contain contributions from older stellar populations. We use simple 
single SSP fits to the SHELS spectra to study the influence of emission lines on the k- 
correction. The effects of emission lines can be ignored for rest-frame equivalent widths 
(REWs) < 100 A depending on required photometric accuracy. We also provide ana- 
lytic approximations to the k-corrections determined from our model fits for z < 0.7 as 
a function of redshift and D„4000 for ugriz and UBVRI {a of the scatter is typically 
~ 0.10 mag, the root-mean-square typically ~ 0.15 mag). Again, the approximations 
using L'„4000 are as good those based on a suitably chosen observed galaxy color. We 
provide all analytical approximations in an online calculator. 

Subject headings: Galaxies 
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Introduction 



The expansion of tlie Universe shifts the spectrum of galaxies redward with respect to the 
observer. Observing galaxies at different cosmological distances in a fixed bandpass thus leads to a 
different absolute magnitude in the bandpass after correcting for luminosity distance. To compare 
the photometric properties of galaxies at different redshifts, one needs to bring the photometry of 
these galaxies onto a common system. The correction that places the absolute magnitude on a 
common system is the k-correction (after iHubbld 119361 ) . The k-correction is usually defined with 
respect to photometric observations of a galaxy in its rest frame, i.e., z = Accurate k-corrections 
have increasing importance as cosmological observations become more precise. Here, we provide 
k-corrections based on Z)„4000 — a spectroscopic indicator that is a measurement of the strength of 
the 4000 A, or Balmer break — and redshift accurate to 0.1 mag with negligible systematic errors. 

The k-correction i^g/j is the correction in magnitudes between the observed magnitude in 
the fi lter R and true rest-frame absolute magnit ude Mq in the filter Q at distance modulus DM(z) 
(e.g., iHogg et ahlbood : Isianton fc RoweislboOTl ^: 



lUR = Mq + DM(z) + Kqr{z 



(1) 



where DM(2;) = 51og(dL/10pc) — 5 and is the luminosity distance. We consider the case where 
filter Q is equal to filter R. 

To determine the k-correction of a galaxy, one approach is to model the spectral energy dis- 
tribution (SED) of the galaxy over the wavelengths necessary to determine the k-correction. The 
quality of the k-correction is then limited by the adequacy of the models; when using photometric 
SEDs rather than spectra, the quality is limited by the number (and width) of the photometric 
bands used. A second approach is direct determination of the k-correction from the observed spec- 
trum of the galaxy. Using observed spectra for calculating the k-corrections removes the need for 
model assumptions and empirically addresses all of the physics required to determine the SED of 
a galaxy. 



Recently, IChilingarian et al.l (|20ld . hereafter ICMZld ) approximated the k-corrections for nine 
optical and near-infrared filters {ugrizY JHK) with analytical functions of only two parameters: 
redshift and observe d galaxy color. They used the measured redshift to fit a set of SSP models 
based on pegase.2 (jFioc Rocca-Volmerangd 119971 ) to the photometry in the nine bands. From 
the fitted models, they derived k-corrections for the filter set. 

Here, we discuss k-corrections determined directly from spectra from the Smithsonian Hec- 
tospec Lensing Surve y. SHELS is a sp ectroscopic survey covering 4 0° on the sky to a limiting 
magnitude R = 21.0 (jGeller et al.ll2005l . Kurtz et al. 2011, in preparation). Instead of observed 
galaxy color, we use the -0^4000 of the galaxy together with redshift to parametrize the k-correction. 



^ To remove lar^e systemat ic uncertainties from the k-correction, sometimes the choice for a non-rest frame is made 
(e.g., iBlanton fc RoweislboOTi '). 
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L>„4000 is the ratio of average flux red- and blueward of the so-called 4000 A break, or Balmer 
break 

MIOO ^ ,^ 

^"4000 = (2) 

J3850 -/^"^ 

where f\ is the rest-frame flux density of the galaxy and each wavelength is the rest-frame wave- 
length. We indicate the regions where we measure -D„4000 (Fig. [1]) with respect to the Chariot 
& Bruzual (private communication 2010) (hereafter CB07) models at three different redshifts. We 
also show the g and r band throughput curves from the SDSS Web sit^. 

The wavelength range covered by the Hectospec spectra determines the choice of fllters for 
which we can determine the k-corrections. We can determine k-corrections only for bandpasses 
that satisfy the "Goldilocks principle" : if the bandpass is too blue, we can determine the rest-frame 
magnitude, not the observed-frame magnitude; if the bandpass is too red, we can determine the 
observed-frame magnitude, not the rest-frame magnitude. We need both observed- and rest-frame 
magnitudes to determine the k-correction. Two fllters that satisfy this principle are the SDSS g 
and r bands. We derive these k-corrections (and analytic approximations) in Section SI 

We discuss the SHELS spectroscopic data in Section [2l In Section [3l we describe the methods 
used to calculate the k-corrections, the sample selection, and the use of -D„4000 as an indicator 
for galaxy type. We present the k-correctio ns in Sec tion HI We compare our results with our SSP 



flts based on CB07 to the previous work of ICMZlOl in Section [5l We discuss the accuracy of our 



L'„4000 in Section [6l In Section [7] we discuss the impact of emission lines on the k-correction. We 
provide the k-corrections for SDSS bandpasses ugriz and the Johnson-Cousins UBVRI derived 
from our model flts to the spectra in Section [5J We summarize our results in Section [21 

Throughout this article we assume a flat universe with Hn = 71kms~^Mp c~^ Q.^ = 0.27, 



and f^A = 0.73. All quoted magnitudes are on the AB-system (jOke &: Gunnlll983l ). 



2. The data 

Here, we describe the observations of SHELS. We flrst describe the imaging and spectroscopic 
data. We then address the effects of flxed solid-angle aperture sizes on our measurements. We 
conclude with the selection criteria of our sample of spectra from which we determine our k- 
corrections. 



^See|http : //www . sdss . org/dr7/ instrvunents/imager/index .html 



2.1. The observations 



We draw our data from an expanded SHELS catalog ( cf. iFabricant et alj l2005l : iGeller et al 



2005 



Kurtz et alJl2007l : lGeller et al 



201C 



Westra et alJl201C : Woods et al 



2010j, who use a previous 



version of the catalog). The main differences with the prior catalog are: (1) we select fainter 
targets, because the survey now reaches R = 21.0, instead of i? = 20.3; (2) we have more redshifts 
available (15,513 unique redshifts compared with 11,447 previously); and (3) we have updated the 
photometry. We next summarize the galaxy catalog construction. We discuss the details in a 
forthcoming article (Kurtz et al. 2010, in preparation). 

We constructe d the SHELS galaxy catalo g from the i?-band source list for the F2 field of the 



Deep Lens Survey (jWittman et al 



2002 



20061 ) . The imaging of F2 is extremely deep; the imaging 



has a 1 o" surface brightness limit of 28.7mag/n" in R. The effective exposure time with the 
MOSAIC I imager (jMuller et al.lll998l ) on the Kitt Peak National Observatory Mayall 4 m telescope 
is about 14,500 s. All exposures a r e obs erved in less than 0'.'9 seeing between 1999 November and 
2004 November. IWittman et al.l (120061 ) describe the reduction pipehne. The 4.2 F2 field is 
centered at a = 09''19™32!4 and 5 = +30°00'00". We exclude regions around bright stars (~ 5 % of 
the total survey), resulting in an effective area of 4.0 D°. We use surface bri ghtness and magnitude 
in combination with information from the latest SDSS data release (DR7; lAbazajian et al.ll2009l ) 
to separate stars from galaxies. This selection removes some active galactic nuclei (AGNs). The 
final catalog contains 20,116 galaxies with a total i?-band magnitude of < 21.0. 



We acq uired spectra for the galaxies with the Hectospec fiber-fed spectrograph (IFabricant et al 



19981 . 12OO5I ) on the MMT on various nights during the period from 2004 April 13 through 2009 
December 15. The Hectospec observation planning software (jRoll et al.l Il998l ) enables efficient 
acquisition of a magnitude-limited sample. 

The SHELS spectra cover the wavelength range of A = 3700-9150 A with a resolution of ~6 A. 
Exposure times ranged from 0.75 hr to 2hr for the lowest s urface brightness objects in the survey. 
We reduced th e data with the stand ard Hectospec pipeline (IMink et al.ll2007l ) and derived redshifts 



with RVS AO (IKurtz &: Mink! Il998l ) with templates constructed for this purpose (IFabricant et al 
2OO5I ). The 1468 unique pairs of repeat observations imply a mean int ernal error of 56 kms ~^ for 



absorption-line objects and 21kms ^ for emission-line objects (see also IFabricant et al.ll2005l ). 



We photometrically calib rate the spectra by sca ling the flux density obtained by the 1'.'5 diam- 
eter fiber to the total R band. iFabricant et al.l (120081 ) describe the technique used for the Hectospec 
spectra in detail. We can photometrically calibrate the data from this fiber-fed spectrograph be- 
cause it has a particularly stable instrument response. This method assumes that the inner l'.'5 of 
a galaxy is representative for the entire galaxy. See Section 12.21 for a more detailed discussion of 
aperture effects. 



SHELS includes 20,116 galaxies to the limiting apparent magnitude. The integral completeness 
of the redshift survey is 96.5 %, 93.8 %, and 77.3 % to i? = 20.3, 20.6, and 21.0, respectively. FigureE] 
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shows the integral and differential completeness of SHELS as a function of total magnitude. 



2.2. Aperture effects 



The Hectospec fibers subtend I'.'S on the sky, sampling only the inner region of a galaxy. For 
a given g alaxy, the frac t ion of light contained in a fixed angular aperture increases with increasing 



redshift. IXewlev et all (|2005l ) study the influence of fixed aperture size on the determination of 



star formation rate, metallicity, and attenuation in a galaxy. Their conclusion is that when the 
aperture contains more than 20% of the light, the central region is a good representation of the 
entire galaxy. Like I Woods et al.l ([20101), we use the -R-band data to compute the fraction of galaxies 
where at least 20 % of the light is contained in the Hectospec fiber aperture. The majority of our 
galaxies fulfill this requirement (Fig. [3]) . 



Fabricant et al.l (120051 ) find no difference between the D„4000 determined from the SDSS spec- 
tra and the -D„4000 from the Hectospec spectra (their Fig. 12). The median redshift of the SDSS is 
z ~ 0.1: at this redshift, at least 20 % of the light of a galaxy is captured in an SDSS fiber (Fig. [3]). 

Galaxies with less than 20 % of their light captured by the Hectospec fibers are predominantly 
at low redshifts. A small difference in color between the nuclear region and the entire galaxy 
only leads to a small uncertainty in the k-correction, due to the low redshift; at low redshifts, the 
4000 A-break has not yet moved out of the filter (Fig. [T]), resulting only in small k-corrections. 

We also assume that the galaxy color measured in a fixed aperture on the sky is a good 
representation of the galaxy. The color obviously varies with aperture in a galaxy with a bulge and 
an extended disk. To investigate the effect of our simplifying assumption, we compare (5 — r) for 
three different apertures. Figure H] shows the comparison between ((7 — r) derived from the SDSS 
fiber spectra covering the inner 3" of a galaxy and from the Hectospec fiber spectra subtending 
the inner 1'.'5 of a galaxy. The figure also shows a comparison b etween the {q — r) from the SDSS 



2001 



Yasuda et al 



Petro sian magnitudes, which include nearly all the galaxy light (iBlanton et al 
2OOII ). and {g — r) determined from the Hectospec fiber spectra. The larger scatter in the lower 



diagram refiects the large uncertainties in the Petrosian magnitudes. 

We find no significant trend for [g — r) within different apertures as a function of redshift. 
There is a slight offset (median of {g - r)fiber,SDSS - id - r)uecto = -0.02 and {g - r)petro,SDSS - 
{9 ~ '')Hecto = —0.05), smaller than the median uncertainty for each of the individual SDSS g and 
r m agnitudes. This offset o ccurs because the outer parts of a galaxy is bluer than the bulge (see 
also iKauffmann et al.ll2003l ). The effect on the determination of the k-corrections is negligible. 



We conclude that D„4000 and [g — r) determined from a Hectospec spectrum are reasonable 
representations of the 1)^4000 and {g — r) obtained from a spectrum that covers the entire galaxy. 
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2.3. Sample selection 

To derive the best k-corrections, we select a subsample of spectra from the sample of 15,513 
galaxies with unique redshifts. We retain the spectra of galaxies satisfying all of the following 
criteria: 

• The spectrum was observed after 2004 October 30. Until this date , the atmospheric disp ersion 



compensation prisms of Hectospec rotated in the wrong direction ([Fabricant et al.ll2008l ). The 
spectra before this date have improper spectrophotometric calibration. However, the redshifts 
for these galaxies are not affected. 

• The model fit (Section 13. 2|) succeeds and the reduced of the fit is less than 10. This 
restriction mostly removes strong broad lined AGN, very strong emission-line galaxies, and 
spectra with bad sky-subtraction. 

• The median of the signal-to- noise ratio (S/N) per pixel for ten regions free of night-sky 
emission lines exceeds 2. 

• The redshift of the galaxy allows the determination of both the observed- and rest-frame 
magnitude. 

After applying these selection criteria, we have samples of 11,707 and 5993 spectra to calculate the 
k-corrections for g and r, respectively. 



3. Method 

Here, we describe the method for determining the k-corrections from our spectra. We first 
outline the formulae used to derive the magnitudes. We also describe the fitting of SSP models to 
our spectra. Finally, we discuss the advantages of using L'„4000 rather than observed color as an 
identifying characteristic of a galaxy. 



3.1. Synthetic magnitudes 

We use fiux-calibrated spectra (and SSP model fits; Section 13. 2p to calculate the observed- 
and rest-frame magnitudes required for the k-corrections. We compute the effective fiux-density 
Fx from the observed spectrum for a particular bandpass 

^ jf{\)\T{X)d\ 
^ f\T{X)dX ' ^ ' 
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where /(A) is the flux density as a function of wavelength A in A, i.e. the flux-calibrated spectrum 
in units of ergs~^ cm~^ and T(A) is the filter transmission curve per photon with wavelengthjf] 



We convert Fx into the effective flux-density as a function of frequency F,^: 

3.34 X IQ-^^Fx 



F, 



\2 



(4) 



where F,y is in erg s ^ cm ^ Hz ""^ , Fx in erg s ^ cm ^ A ^, 
angstrom. We calculate the artificial magnitude as 



and AcE = / AT(A)dA/ / T{X)dX is in 



m 



-2.5 log 



Fu 



3631Jy 



(5) 



where 3631 Jy is the zeropoint of the AB-magnitude system (|Qke Gunnlll983l ). 

Throughout this article, unless otherwise noted, we determine the k-corrections from synthetic 
observer and rest-frame magnitudes. We obtain the observed colors of the galaxies from the SDSS 
fiber magnitudes; when these magnitudes are unavailable, we determine them from our synthetic 
magnitudes. 



3.2. Model fitting 



We use the method of lTremonti et al.l (|2004l ) to fit a linear combination of 10 SSP models to 
each spectrum. The SSP models are from CB07. These models include an im proved treatment of 
the th ermal pulse asympt otic giant branch phase with respect to the models of lBruzual &: Chariot 
(|2003l ) (see iBruzuall l2006l . for more details). The 10 different ages for the populations are 0.005, 
0.025, 0.1, 0.25, 0.5, 1, 1.4, 2.5, 5, and lOGyr. The models have solar metallicity. Figure [T] shows 
the indi vidual model spectra normalized to 1 at 5500 A toget her with the SPSS g a nd r bands. We 



use the lCalzetti et al.l (2000ll extinction law rather than the ICharlot &: Falll (j2000l ) extinction law 
used bv lTremonti et al.l (|2004l ). but the difference is small (10 % at the extremes of the wavelength 



range 4000 A - 10,000 A). 

We fit the models to the wavelength region between 4000 and 8500 A; we use the fits to deter- 
mine L'n4000 (Section 13. 3p and remove residuals of night-sky features in our spectra (Section H]). 
The cut at the blue end of the spectrum is due to the uncertain calibration of Hectospec in that re- 
gion. Beyond 8500 A light emitted from the home sensors in the fiber positioner of the spectrograph 
affects the spectrum. Because this effect is additive, we cannot remove this problem. 



Blanton fc Roweia (|2007l ) note that many authors give the transmission curve as the contribution to the detector 



signal per unit of energy (T'(A)), rather than per photon (r(A)). They relate to each other as r(A) oc T'{\)/X. 
However, accidentally using T'{X) as r(A) has a very small impact (at the very worst 0.07 mag, and typically 0.02 mag 
or less for this survey). 
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Figure [5] shows the predicted k-corrections as a function of redshift for each model without 
attenuation. We include two models with attenuation to show its effect on the k-correction. We 
apply a typical reddening for star-forming galaxies to the model with the youngest stellar population 
(Ay = 1) and a reasonably extreme reddening for early- type galaxies to the model with the oldest 
stellar population {Ay = 0.5). Figure [5] shows approximately the full range of expected k-corrections 
modulo the effects of emission lines (see Section [7]) . 



3.3. i:>„4000 versus color 



We use -0^4000 as a meas ure of galaxy type because it is a n indicator of the a ge of the 
stellar population (see Table [U iBalogh et al.l Il999l : iBruzuall Il983l ) . IWoods et al.l (j2010l ) examine 
the fraction of galaxies with emission lines as a function of 1)^4000 (they remove AGNs from their 
sample). Galaxies with a low L'„4000 invariably show e mission lines; the f ractio n of galaxies with 
emission lines decreases rapidly with increasing Z)n4000. iKauffmann et al.l (120031) show a very clear 



Kauffmann et al. also show 



relation between Z)„4000 and (a measure for the REW of H5). 
clear evolution in both Z)„4000 and H(5^ after an instantaneous burst of star formation as a function 
of time. 

Figure |6] shows the distribution of L'„4000 for our sample, along with the fraction of galaxies 
with at least one emission line. We classify a galaxy as having an emission line when the absorption- 
corrected REW is at least 5 A for any of the following emission lines: Ha, H/3, [O iii] A4959, 
[O III] A5007, or [On] AA3726,3728. Our distribution is double-peaked w ith a local minimum around 
L'„4000 = 1.46 {thick dashed line), very close to the lWoods et al.l (|2010l ) value of 1.44 (dotted line). 



Almost all galaxies with a low L'„4000 (< 1.2) show emission lines indicating an episode of 
active (or very recent, roughly less than lOMyr ago) star formation, yielding a spectrum dominated 
by young, hot stars. The fraction of galaxies with emission lines gradually decreases with increasing 
-D„4000. For galaxies with a large 4000 A-break, an old stellar population dominates the spectrum. 
In our sample, toward the highest Dri,4000 the number of galaxies with emission lines increases 
slightly because of the presence of AGNs. 

An advantage of Z)„4000 over observed color is that L'„4000 is insensitive to reddening; L'„4000 
is measured over a small wavelength range (see Fig. [1]). In addition to the Dn^OOO for the 10 SSP, 
we show the L'„4000 for two of these models with some attenuation in Table [TJ Attenuation affects 
L'„4000 only minimally. Furthermore, the small wavelength range needed to calculate Dn'iOOO 
allows the use of uncalibrated spectra, provided that the sensitivity of the spectrograph does not 
vary significantly over the relevant wavelengths. 

Another advantage of L'rt4000 is — by its definition — independence of redshift; the observed 
color of a galaxy changes dramatically with redshift. At different redshifts, the bandpasses probe 
different parts of the spectrum. Thus, the k-corrections are not (necessarily) equal for different 
filters. Consider, for example, the 4000 A break, a large-scale feature (in wavelength) in the galaxy 
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spectrum where the flux can change significantly. The 4000 A break is just in the g band at z ~ 0, 
leaves it around z ~ 0.35 and enters the r band, and leaves the r band at z ~ 0.75 (see Fig. [T])jl| 

Figure [7| shows observed {g — r) as function of -D„4000 for both the separate SSP models and 
the individual galaxies. We indicate the median color of the galaxies binned by redshift as function 
of Dn4000. Toward higher D„4000, {g — r) becomes less sensitive than L)„4000 to the age of the 
stellar population of a galaxy. Furthermore, the models show that {g — r) is not a monotonic 
function of redshift; L)„,4000 is redshift-independent. 

A spectrum of a galaxy adequate to determine its redshift (a necessary quantity to calculate 
the k-correction) usually straightforwardly provides the L>„,4000 of that galaxy. The -D„4000 is 
more robust for galaxies where older stellar populations dominate the spectrum; the redshift is 
then usually determined from the Ca H+K lines. 

We use the -D„4000 determined from our model fits. This approach allows us to determine the 
1)^4000 more accurately for lower S/N spectra, because we use almost the entire wavelength range 
of the spectrum to constrain the model fit; the fits refiect the noise over the entire spectrum rather 
than the noise in a small portion of the spectrum around 4000 A. 



4. Empirically determined k-corrections 

Here, we derive the empirically determined k-corrections from our spectra and from the model 
fits to the spectra as a function of redshift and 1)^4000. We also provide an online calculatorlf] 

To minimize the impact of large residuals of night-sky features in our spectra, we replace the 
observed spectrum around these wavelengths with the model fit. We follow the same procedure 
for observed wavelengths blueward of 4000 A and redward of 8500 A because of the uncertainties 
in the flux calibration of our spectra. Using the deflnition of the k-correction (eq. [1]) and the 
prescription for determining artiflcial magnitudes from our spectra (Section EH]), we determine the 
k-corrections from both our spectra and the model flts. Figure [8] gives the resulting k-corrections 
for the g and r band. Each galaxy is color-coded by its -D„4000. 

To use these results for other observed galaxi es, we ap proximate these results with a two- 



dimensional surface in a fashion similar to that of ICMZld . We use a minimization to fit a 
two-dimensional surface to our data. We describe the fitting method in detail in the Appendix. 
We give a short summary here. To derive the fit, we use the redshift, -D„4000, and determined 
k-correction (respectively, Xk, Uk, and Zk in eq. |A2j and eq. |A3j ) for all galaxies satisfying the 
sample selection of Section [2.31 to derive the fit. We choose the maximum polynomial orders for 



^Similarly, the Lyman break causes drastic changes in galaxy colors as a function of redshi ft, which is a property 
employed in selecting higher redshift galaxies with the drop-out technique I Steidel et al. 1995l l. 



^See |http : //tdc-www. cf a. harvard. edu/iiistruments/hectospec/progs/EOK/[ 
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the redshift and Z?„4000 {Nx and Ny in eq. |Alj . respectively) to be the lowest orders that do 
not show an obvious pattern in the residuals after fitting. If we use too many orders, the fit 
provides inaccurate results for the highest redshifts or extreme L'n4000 values, because the number 
of galaxies constraining the fit is limited. At 2; = the fit is constrained by the definition that the 
k-correction is zero (Section [T|). For our fits, = 3 and Ny = 3 are appropriate. 

Tables with the coefficients for the approximation are in the Appendix (Tables IA.mA.4]) . Fig- 
ure [9] shows the residuals as a function of L'„4000 and redshift for g and r. The thick colored lines 
indicate the median of the residuals binned by Dn4000 to show the absence of structure in the 
residuals as a function of redshift. Table [2] contains statistics {a of the Gaussian fit, rms, and the 
68.3 % range around the mean) on the residuals for the fit. 



We also determine the coefficients for our fit as a functi on of (q — r) (now yk in eqs. |A2 



and jA3j ) and redshift in order to compare our results with ICMZld (see Section [5]) . We show 
the residuals in Figure [TOl for g and r. The coefficients for the approximation are in the Appendix 
(Tables [A.SIIA.Sp . Table [2] contains statistics on the residuals. We note that the analytic approx- 
imations using on D„4000 show the same, or lower, residuals for the r and g bands, respectively, 
compared with those approximations of the same polynomial order based on (g — r). 



5. Comparison with previous work 



The advantage of using k-corrections determined directly from spectra relative to those from 
models or model fits is that we make no assumptions about galaxy evolution, attenuation laws, 
metallicities, orientation, etc. Using the spectra directly also eliminates the need for constructing a 
set of models that — in linear combination — represent fundamental galaxy types. However, models 
do not have the (photon)noise that the observations contain. Here, we compare our empirical 
k-corrections with the models of CMZIC . 



Figure [8] shows the model tracks of pure SSP models from Figure [5] plotted over our data. The 
range of spectroscopically determined k-corrections matches the range of the models at the high 
-0^4000 end extremely well. At the low L'„4000 end, there are few galaxies with k-corrections as 
low as predicted by the youngest SSP models without attenuation. This difference implies that 
these low D„4000 galaxies are always (significantly) attenuated and/or these galaxies rarely have 
a stellar population that consists purely of young stars. 



We next compare our k-corrections with those from the literature. ICMZIG provides k-corrections 
as a function of redshift and just one other parameter, the observed galaxy color. They use pho- 
tometry in up to nine optical and near-infrared bands for galaxies with known redshift to fit a 
grid of SSP models from pegase.2 (jFioc fc Rocca-Volmerangd 119971 ). In Figure [TT] we show the 
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k-corrections for our galaxies determined with the prescription of CMZlol l^ 



Fig ure 1121 shows the difference between the k-corrections derived from our data and those from 
CMZld for SHELS galaxies as a function of both (g — r) and D„4000. The k-corrections for the r 
band agree extremely welL The k-corrections for the g band also agree well up to z < 0.3 (68.3 % 
range around the median difference is 0.15 and 0.12 for g and r, respectively at z = 0.3). For larger 
z, they start to disagree for galaxies with large 1)^4000. The cause for this difference is twofold. 



At a redshift of ~ 0.3, the iCMZld sample does not contain as r nany green-to-blue gal axies 
as red galaxies. They drew their sample of galaxies from SPSS DR7 (lAbazajian et al.ll2009l ) and 
UKIRT Infrared Deep Sky Survey DR5 (iLawrence et al.l 120071 ) with redshifts between 0.03 and 
0.6. The galaxies in SDSS DR7 have spectroscopy for at least two distinct samples: the main 
galaxy sample (r < 17.77) with a median redshift around 0.10 and the luminous red galaxy sample 
(r < 19.5) selected with the majority of the galaxies at redshifts higher than 0.15. 

In contrast, SHELS reaches much fainter than SDSS (i? < 21) and selects galaxies based only 
on their total magnitude. Thus, SHELS is not biased toward or against a particular type of galaxy 
(except against AGN due to their stellarlike mor phology) . SHELS contains more (in number) 
intrinsically bluer galaxies at higher redshifts than lCMZlC . 



CMZld only selects galaxies with redshifts up to 0.6. Thus, their fit is only valid to approxi- 
mately that redshift; SHELS contains galaxies to z = 0.7. 



However, the differences with lCMZlG for the galaxies with the largest -0^4000 require further 
discussion. Figure [13] shows the differences between the magnitudes derived from the spectra and 
the models in the observed and rest frames. Modulo the presence of emission lines (Section [7| 
and noise, these magnitudes should be the same. Apart from the panel with observed gf-band 
magnitudes, all other panels show that the magnitudes derived from the spectra and the models 
are consistent with each other. 

We investigate the overestimation in flux indicated by Figure [13] of the high- Z)„ 4000 galaxies 
from the models with respect to the spectra by combining the spectra of these galaxies. We note 
that for the galaxies with the oldest stellar populations there is less flexibility in fitting the models 
than for galaxies with younger populations because only old SSP models can fit; for younger popula- 
tions conceivably different combinations of young and old SSP models (and an increased/decreased 
amount of attenuation) can yield a good fit. 

Figure [14] shows the averaged luminosity-weighted rest-frame spectrum — of both the actual 
spectra and the model fits — of 84 galaxies with 0.4 < z < 0.6, L'n4000 > 1.7, and a difference 



^ Blanton fc Roweis (2007) also p rovide a method to determine the k-corrections. The method of CMZlol is similar 
to that of iBlanton &: Roweid (|2007l ) in the sense t hat both use template sp ectra to fit the photome tric data, but 
the template set of CMZlol differs from the set of Blanton &: Roweii (2007). Furthermore, CMZ1(]| do a detailed 
comparison and conclude that the results from both approaches are similar. 
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between the magnitude derived from the spectrum and the fit that is larger than 0.2. We also 
show the difference between the averaged spectrum and averaged model in the middle panel and 
the relative difference between the averaged spectrum and averaged model with respect to the 
averaged model. 

The models deviate from the spectra at the bluest rest-frame wavelengths. The relative differ- 
ence there is large. Because the features in the average fit also appear in the spectrum, the flux of 
the average spectrum goes only marginally below 0. These rest-frame wavelengths correspond to 
the bluest part of the spectrum (where the spectrograph is not very sensitive) and it appears that 
the sky is somewhat oversubtracted. We suggest the use of the k-corrections determined from the 
model fits for the galaxies with high L'n^OOO at the highest redshifts. 

Overall, we fi nd very g ood agreement between our k-corrections us ing the full spectrum and 



the prescription of ICMZIQ based on SSP fits to photometry. Similar to lCMZlQ (who use redshift 
and observed color), we can approximate the k-correction with only two parameters: redshift and 
-D„4000. There are no systematics remaining in the residuals of our approximation. 



6. Uncertainties 

In Section [3.31 we determined -0^4000 from the model fits to the spectra. Not every spectrum 
is suitable for fitting SSP models; in those cases one can determine 0^4000 from the spectrum. 
In this section we consider the effects of using the spectrum-derived -D„4000 on the determination 
of the k-correction instead of using the model-derived -D„4000. We compare the accuracy of the 
model-derived -D„4000 with the spectrum-derived -D„4000. 

The accuracy of L>„4000 is obviously a function of the S/N of the spectrum. The S/N for each 
of our spectra is the median of the S/N per Hectospec pixel measured in 10 wavelength ranges in 
the observed spectrum that are free of night-sky emission and sample the range over which we fit 
our models to the spectrum (4100-4300, 5000-5200, 5470-5570, 5700-5850, 6000-6200, 6400-6600, 
6650-6700, 7000-7200, 7600-7700, and 8050-8250 A). Figure [E] shows the relative difference of the 
spectrum- and model-derived D„4000 as a function of S/N. The accuracy of Dn^QOO increases, as 
one might expect, with increasing S/N. 

Figure [15] also shows the distribution of the relative difference between the model- and spectrum- 
derived -D„4000 binned by S/N. The 68.3% range around the median increases with decreasing 
S/N (0.21, 0.13, 0.11, and 0.07, respectively, for the four S/N ranges indicated). The values of the 
higher S/N spe c tra ar e consistent with the comparison between SDSS and Hectospec spectra of 



Fabricant et al.l (|2008l . rms = 0.09). 



Figure [15] also shows the distribution of model-derived Z)„4000 (because we use this throughout 
this article) in the same S/N bins. Galaxies with the most poorly determined -D„4000 are typically 
the low surface brightness — also intrinsically bluer — galaxies. 
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Table [3] shows the k-corrections and their uncertainties for a few example galaxies given the 
uncertainties in Z)„4000 and {g — r). The three uncertainties for Z)„4000 correspond to the 68.3% 
range around the median for very low, moderate, and high S/N spectra. The uncertainties in the 
corresponding color are those in the SDSS Petrosian and fiber colors. The column at indicates 
the corresponding range in the calculated k-correction. We choose the redshifts of the galaxies, 
z = 0.30 and 0.55, close to the maximum redshifts where we can determine k-corrections for the 
r and g bands, respectively. We choose 1)^4000 values that are typical for the two populations of 
galaxies (Fig. L»„4000 = 1.25 and 1.65). 

The D„4000 from a spectrum with S/N > 2 (cj£)„4ooo ^ 0.11) yields an uncertainty in the k- 
correction, similar to the uncertainty in the k-correction determined from {g — r) given the typical 
uncertainties in the galaxy color. For the Dri4000 from a spectrum with S/N < 2 {(TD„iooo > 0.11), 
the uncertainty in the k-correction is somewhat larger than the uncertainty in the k-correction 
determined from (g—r). In most of the examples, the uncertainty in the SDSS magnitude dominates 
the uncertainty in the final absolute magnitude, not the uncertainty in the k-correction. 



7. Emission lines 

Emission lines with large REW can significantly affect the flux density measured through a 
filter0 Figure [13] shows the influence of the emission lines on the magnitudes. The galaxies with 
low Z)„4000 and |mspec — f^modl ^ are almost all emission-line galaxies. 



As IWoods et al.l (120101 ) show (see also Section 13.31 and Figure [6]) , nearly all galaxies with low 
Z?„4000 have emission lines. These galaxies make up a large fraction of our sample. We thus 
investigate the influence of the emission lines on the derived k-corrections. 

We examine five galaxies with differing emission-line REWs. Figure [16] shows the spectra of 
these galaxies. We calculate the magnitudes and k-corrections as a function of redshift from the 
rest-frame spectra and the model fits of these galaxies. We assume that the difference between the 
spectrum and the model results only from the emission lines. We plot the difference in derived 
magnitudes and k-corrections for both bands in Figure [TtI^ 

The entrance and departure of emission lines from the filter as a function of increasing 
redshift dominate the curve that traces the difference in magnitudes. At z = the H/3 and 
[O III] AA4959,5007 complex are in the g filter and leave it with increasing redshift. Around z ~ 0.1 
the [On] AA3726,3728 doublet enters the filter, causing the difference in the magnitude to increase 



^Narrowband surveys employ exactly this feature to increase the contrast between observations done with a very 
narrow filter and a broadband fil ter; the narrower th e filter becomes, the larger the difference in magnitude between 
the two filters (e.g., Waller f 99ol: jPascual et al.llioOTi '). 

*The spectra of these galaxies normally evolve over time. The point of this figure, however, is to demonstrate the 
infiuence of emission lines, not galaxy evolution. 
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only slightly. For the r band, the rippling shape of the curve results from [S ii] AA6733,6718, 
[N il] AA6550,6585, and (predominantly) Ha leaving the filter with increasing redshift and having 
left around z ~ 0.1, while the H/? and [O iii] complex starts to enter the filter at the same time. 
Around z ~ 0.35 the complex starts to leave, while H7 enters. At z ~ 0.55 H7 leaves and H(5 
enters. 

FigurefTTlshows that (1) the influence of emission lines is clearly correlated with their equivalent 
width, and (2) apart from galaxies with extremely large REW emission lines (REW ^ lOOA), the in- 
fluence of emission lines on the broadband g and r magnitudes — and, consequently, k-corrections — is 
reasonably small. The significance depends on the accuracy of the photometry and/or the required 
precision of the k-corrections. F or an accuracy < 0.1 mag k-corrections need to take emission-line 



flux into account. For example, iBrown et al.l ([2003) use the influence of strong emission lines on 



the observed galaxy colors to select extremely metal poor galaxies. 



8. Additional k-corrections 

In Section Owe determined the k-corrections for the SDSS g and r bands; here we consider k- 
corrections in other photometric bands. Because emission lines are important for only the strongest 
emission lines in a galaxy spectrum and our model fit extends beyond the wavelength range of the 
spectrum, we can calculate the k-corrections for any optical bandpass over the redshift range where 
we sample enough galaxies: i.e., z < 0.7. 

Here, we derive the k-corrections as a function of redshift and .Dn^OOO fo r the SDSS bandpasses 



ugriz and the Johnson-Cousins UBVRI given in Table 2 of Bessell (jl99ol ). Figure [18] shows the 
10 normalized bandpasses we used to calculate the k-corrections. Figure [T9l showcases the results 
for the SDSS u and the Johnson-Cousins V band. The figures for these and the remaining filters 
are also available onlinelfl Tables IA.9tiA.18l contain the coefficients for the analytic approximation 
using redshift and Z)„4000, Tables IA.19IIA33] contain the coefficients for the approximations using 
redshift and color (SDSS filters only; we do not have imaging available in the Johnson-Cousins 
bandpasses). We give the statistics of the residuals in Table [4] to indicate the quality of the fits. 



For the ugriz bandpasses, our k-corrections match th ose oflCMZld at z < 0.6 for galaxies with 



high Z'„4000 and at z < 0.3 for galaxies with low L'„4000. ICMZlOl selects galaxies with redshifts up 



to 0.6 and their number of green galaxies decreases steeply beyond red shift 0.3 . Furthermore, toward 



redder filters the differences between our k-corrections and those of lCMZld decrease; because the 
4000 A break goes through the redder filters at higher redshifts, the k-corrections for the redder 
filters are also smaller than for the bluest filters. 

The residuals from the analytic approximation using Z)„4000 and redshift are about the same 



^See |http: //tdc-www. cf a. harvard. edu/ instrument s/hectospec/progs/EOK/figures . shtml[ 
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as or smaller than those using observed galaxy co lor and r edshift. This result is, however, not sur- 
prising; the optical observed galaxy colors used bv lCMZld coarselv probe Z)„4000 for the particular 
filter where they calculate the k-correction. 



9. Summary 

We use the 15,513 SHELS spectra to determine the k-corrections for the SDSS g and r bands 
empirically. These k-corrections have the major advantage that they are model-independent. We 
express the analytic approximations for the k-correction as a functio n of two parameters: redshift 



and L)„4000. These parameters have some advantages over those from lCMZlQ : redshift and galaxy 
color. The advantages include (1) -D„4000 is redshift-independent, (2) L>„4000 is an indicator of 
the age of the stellar population of the galaxy, (3) D„4000 barely suffers from attenuation, (4) the 
spectrum used to determine the redshift of a galaxy can also yield Z)ri4000, and (5) the spectrum 
does not necessarily have to be flux-calibrated to determine 1)^4000. Furthermore, the analytic 
approximations based on D„4000 have about the same (r band) or less scatter {g band), compared 
with the approximations using (g — r). The a of the scatter in the residuals is ~ 0.08 mag, the rms 
is ~ 0.15; the residuals also show no systematic bias. 

From CB07 SSP models we indicate the range of k-corrections that we determine from observed 
galaxies. The models with the youngest and oldest stellar populations with a small amount of 
attenuation span the entire range of determined k-corrections. We conclude that galaxies are 
always (significantly) attenuated and/or rarely have a stellar population that consists purely of 
young stars. 



We compare our results to the analytic approximations given bv ICMZld . Up to z < 0.3 our 



results agree well for both bands. Differences between the two surveys arise due to the fa ct that 



SHELS goes to significantly fainter magnitudes than SDSS (the survey from which ICMZlOl draws). 
Furthermore, our total magnitude selection criterion allows us to sample the entire range of galaxies 
up to z ~ 0.7, rather than selecting for luminous red galaxies (like SDSS). We make the analytic 
approximations available in an online calculator (see footnote 7). 

We check the influence of emission lines on the k-corrections. The influence of the emission 
lines depends on the REW of the line; the larger the REW of the line, the larger the influence 
on the k-correction is. Depending on the accuracy of the k-corrections required, the influence of 
emission lines can typically be ignored, except those few cases in which the REW is very large 
(REW > 100 A gives Am > 0.10 mag). 

Finally, we provide k-corrections as a function of redshift and Z)„4000 using our model fits 
for SDSS ugriz and Johnson-Cousins UBVRI. We determine these k-corrections over the redshift 
range where we sample enough g alaxies: i.e., z < 0.7. We compare our results to the analytical 



approximations given by ICMZlOl and sh ow that w e cover a larger redshift range, where we sample 



more (in number) "green" galaxies than ICMZlOl do. For spectra with S/N > 2 {(Td„4:000 ^ 0.11 
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the k-corrections based on D„4000 are as good as those based on color. 

Despite the fact that the scatter in our k-corrections is typically ~ 0.1 mag, the uncertainties 
in the observed magnitude are of a similar size or larger. Thus, uncertainties in the magnitudes 
dominate the uncertainty in the final, k-corrected magnitude. 

This empirical approach to determine k-corrections can be extended to the near- and mid- 
infrared to include spectroscopic features, such as the Paschen and Brackett hydrogen lines (and 
breaks) and polycyclic aromatic hydrocarbon features. The technique can also be extended to 
higher redshifts where the 4000 A break moves through the reddest optical filters and into the 
near-infrared. 
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In Section [5] we describe the empirical k-correction as a result of redshift and either D„4000 
or color. Here, we describe the method used for the fitting and give the coefficients for the approx- 
imations. 

We use minimization to determine the coefficients for our analytic approximations. We 
describe our method in the context of fitting a two-dimensional surface. The surface z as a function 
of coordinates x and y is given by 



Acknowledgments 



A. K-correction approximations: surface fitting 
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where Cij are the coefficients and and Ny are the maximum polynomial terms for x and y, 
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respectively. 

The for this surface for Nq data points is defined as 
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The minimization d'x^ /datj = — 2^^(zfc — z{xk,yk))x\yl. = for each coefficient aij yields 



Zkxlvi = z{xk,yk)xlyl. 



(A3) 



We solve these equations using lower/upper (LU) decomposition. For LU decomposition we 
need to write all equations from eq. jA3j as a linear system Ax = B, where x is a vector containing 
the coefficients. Let Sij = "^k^kyi ^-i-j = Ylik^kX^y^, then the linear system we solve looks 
Uke 
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We use the IDL routine la_linear_equation based on LAPACK (jAnderson et al.lll999l ) to solve 
for X. 

We apply this method to the data presented in this article. The variables x^, Vk, and Zk 
represent the redshift, -D„4000 (or {g — r), depending on which of the two we fit), and the k- 
correction, respectively, for each individual galaxy. Ng is the number of galaxies. We show the 
residuals of our fits in Figures [9l [TOl and [19] and online (see footnote 11). 



Tables [X]TVlA. 23] contain the coefficients for the fits described in this article (similar to lCMZlOl ). 
One calculates the k-correction by taking the index of the row multiplied with the index of the 
column (i.e., z*(-D„4000)-', or z^{g — ry), and multiplying that with the coefficient (aij) in the table; 
one needs to do this for all coefficients in the table. The sum of those values is the k-correction for 
the given redshift and L'„4000 (or {g — r) where appropriate). Note that no table contains a row 
with z^ because we demand, by definition, that at redshift z = the k-correction is (Section [1]). 
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Table 1. -D„4000 determined from the models. 



Age 


Av 


L>„4000 


5Myr 


0.0 


0.94 


25Myr 


0.0 


0.97 


100 Myr 


0.0 


1.05 


250 Myr 


0.0 


1.13 


500 Myr 


0.0 


1.26 


1.0 Gyr 


0.0 


1.41 


1.4 Gyr 


0.0 


1.45 


2.5 Gyr 


0.0 


1.58 


5.0 Gyr 


0.0 


1.77 


10.0 Gyr 


0.0 


1.96 


5 Myr 


1.0 


0.98 


10.0 Gyr 


0.5 


2.00 



Note. — These models 
are indicated in Fig. [5l 
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Table 2. Statistics for the analytic approximations for g and r. 









(9- 


r) 


filter 


a 


rms 


68.3 % a 


rms 


68.3 % 








Spectra 






9 


0.08 


0.17 


0.11 0.12 


0.25 


0.17 


r 


0.08 


0.12 


0.09 0.07 


0.11 


0.08 








Models 






9 


0.09 


0.14 


0.11 0.12 


0.21 


0.16 


r 


0.07 


0.10 


0.08 0.06 


0.09 


0.07 


Note. 




a is the 


a of the Gaussian 


fit to the 



residuals for the analytic approximations to the k- 
corrections derived from the spectra and the models 
for g and r with z < 0.68 and 0.33, respectively, as 
a function of Dn4000 and {g — r); rms is the rms of 
the residuals; 68.3 % is the range around the mean 
of the residuals. Figures [9l and [TOl show the residuals 
and the Gaussian fits to the residuals. 



Table 3. Example uncertainties for typical values and uncertainties in D„4000 and corresponding [g — r). 



filter Z Dn^OOO Cpetro CTfib konWOO CrDn4000 CTfc (TDniOOO CTk CTnniOOO CTfc (ff ~ ^) ^(g-r) Cr(B-r),petro CTfc f(g-r),flb ffc 





0.30 


1.25 


0.15 


0.09 


0.61 


0.21 


+0.33 
-0.61 


0.11 


+0.20 
-0.27 


0.07 


+0.13 
-0.16 


0.97 


0.61 


0.17 


+0.16 
-0.16 


0.10 


+0.10 
-0.09 




0.30 


1.65 


0.18 


0.11 


1.15 


0.21 


+0.20 
-0.18 


0.11 


+0.09 
-0.09 


0.07 


+0.06 
-0.06 


1.48 


1.12 


0.20 


+0.10 
-0.15 


0.11 


+0.06 
-0.08 


r 


0.30 


1.25 


0.09 


0.06 


0.21 


0.21 


+0.14 
-0.33 


0.11 


+0.09 
-0.14 


0.07 


+0.06 
-0.08 


0.97 


0.21 


0.17 


+ 0.08 
-0.10 


0.10 


+0.05 
-0.06 


r 


0.30 


1.65 


0.08 


0.04 


0.41 


0.21 


+ 0.03 
-0.04 


0.11 


+0.01 
-0.01 


0.07 


+0.01 
-0.01 


1.48 


0.42 


0.20 


+0.02 
-0.05 


0.11 


+ 0.02 
-0.02 




0.55 


1.25 


0.29 


0.15 


1.08 


0.21 


+ 0.65 
-0.99 


0.11 


+0.37 
-0.47 


0.07 


+0.25 
-0.28 


1.04 


1.29 


0.32 


+0.45 
-0.61 


0.17 


+0.25 
-0.27 


5 


0.55 


1.65 


0.47 


0.24 


2.24 


0.21 


+0.41 
-0.44 


0.11 


+ 0.21 
-0.22 


0.07 


+0.13 
-0.14 


1.65 


2.06 


0.49 


+0.24 
-0.69 


0.25 


+0.18 
-0.27 



Note. — cTpetro and aah are the typical for the SDSS Petrosian and fiber magnitude of a galaxy with chosen L'„4000 and redshift z. The uncertainties aD„4ooo are 
those typical for very low, moderate and high S/N spectra. (Jk is the range in the calculated k-correction corresponding to the uncertainty in Z)„4000 (or {g — r)). 
{g — r) is the typical galaxy color for the chosen i3„4000 and redshift, and cr(g_^)_petro and cr(^g-r),fih are the typical uncertainties in the SDSS Petrosian and fiber color, 
respectively, for such a galaxy. 
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Table 4. Statistics for tlie analytic approximations for UBVRI and ugriz. 







i:»n40oo 






Color 




filter 


a 


rms 


68.3 % 


color 


a 


rms 


68.3 % 


u 


0.16 


0.30 


0.22 


[u — r) 


0.30 


0.48 


0.55 


9 


0.09 


0.11 


0.14 


is -r) 


0.12 


0.22 


0.16 


r 


0.10 


0.15 


0.12 


is -r) 


0.09 


0.17 


0.12 


i 


0.10 


0.15 


0.12 


{9 - i) 


0.09 


0.14 


0.11 


z 


0.09 


0.14 


0.11 


(r — z) 


0.09 


0.14 


0.11 


U 


0.16 


0.30 


0.22 










B 


0.10 


0.17 


0.13 










V 


0.10 


0.14 


0.12 










R 


0.10 


0.15 


0.12 










I 


0.09 


0.15 


0.12 











Note. — The table headings are 

the same as Table [2j The figures at 

|http: //tdc-www. cf a. harvard. edu/ instrument s/hectospec/progs/EOK/figures . shtml 

show the residuals and the Gaussian fits to the residuals. 
Figure [19] shows the figures for u and V . 
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Table A.l. Coefficients for tfie SDSS g band as a function of redshift and Z)„4000 with 

k-corrections derived from the spectra. 



(£»„4000)0 


(£>„4000)i 


(Dn4000)2 


(L»n4000)3 


29.9754 


-78.4292 


64.0436 


-15.8884 


-368.572 


818.057 


-581.89 


134.427 


490.834 


-1084.35 


768.318 


-176.591 



Table A.2. Coefficients for the SDSS g band as a function of redshift and £)„4000 with 

k-corroctions derived from the models. 



(Z)n4000)0 


(Z)n4000)i 


(i^„4000)2 


(L»„4000)3 


6.6812 


-27.6858 


27.7299 


-7.319 


-231.313 


517.183 


-364.586 


82.7252 


330.803 


-735.254 


517.408 


-117.366 



Table A. 3. Coefficients for the SDSS r band as a function of redshift and Z)„4000 with 

k-corrections derived from the spectra. 







(i:'„-iooo)-^ 


(D,, 1000)2 


{D„ 1000)3 


z^ 


25.4725 


-52.2513 


35.6903 


-7.85585 


z2 


-177.822 


336.276 


-203.045 


39.6419 


z^ 


66.6927 


-82.7486 


7.78744 


9.65736 



Table A.4. Coefficients for the SDSS r band as a function of redshift and D„4000 with 

k-corrections derived from the models. 



(Z)n4000)0 


(Z)n4000)i 




(Z)„4000)3 


z^ -2.4549 


3.08855 


-0.85258 


0.169083 


31.7806 


-96.6607 


92.3012 


-26.9024 


z^ -334.912 


776.204 


-593.901 


148.071 
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Table A. 5. Coefficients for tfie SDSS g band as a function of redshift and {g — r) with 

k-corrections derived from the spectra. 









(g rf 


ill rf 




0.126452 


0.0619892 


6.90796 


-2.98031 


z2 


3.51286 


-9.40049 


-9.17117 


6.52303 




-8.56012 


24.3121 


-2.88499 


-3.57257 



Table A. 6. Coefficients for the SDSS g band as a function of redshift and {g — r) with 

k-corroctions derived from the model fits. 





{g-rf 


{g - r f 


{g - rf 


{g - rf 




-0.160448 


1.12881 


5.93095 


-2.68308 




4.04345 


-11.8962 


-6.0987 


5.37315 


^3 


-9.01521 


26.6088 


-6.50937 


-2.17746 



Table A. 7. Coefficients for the SDSS r band as a function of redshift and {g — r) with 

k-corrections derived from the spectra. 





io rf 


{g-rf 


{g rf 


{(] - rf 


z^ 


2.30349 


-7.69531 


11.1784 


-4.33398 


z^ 


-20.1966 


77.9084 


-88.5243 


31.1532 


z^ 


21.5438 


-130.494 


156.555 


-54.8953 



Table A. 8. Coefficients for the SDSS r band as a function of redshift and {g — r) with 

k-corrections derived from the model fits. 





[g-rf 


{g - rf 


{g - rf 


{g - rf 


z^ 


-0.823562 


1.8179 


1.8466 


-1.21972 


z^ 


6.54473 


-13.4719 


3.54192 


1.11608 


z^ 


-28.7761 


59.0299 


-39.6835 


8.96288 
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Table A.9. Coefficients for the SDSS u band as a function of redshift (< 0.7) and D„4000 with 

k-corrections derived from the models. 



(L»„4000)0 


(L>„4000)i 


(L»„4000)2 


(L»„4000)3 


2.43154 


-16.3322 


20.0889 


-6.24137 


^2 -52.0421 


120.631 


-101.738 


31.5506 


147.367 


-326.482 


247.473 


-66.0534 



Table A. 10. Coefficients for the SDSS g band as a function of redshift (< 0.7) and D„4000 with 

k-corrcctions derived from the models. 



(L>„4000)0 


(Z)n4000)i 




(L'„4000)3 


z^ 7.12187 


-28.8064 


28.6961 


-7.57514 


^2 -234.056 


523.817 


-370.133 


84.1717 


333.917 


-742.667 


523.643 


-118.999 



Table A.ll. Coefficients for the SDSS r band as a function of redshift (< 0.7) and Dn^OOO with 

k-corrections derived from the models. 



(L»„4000)0 


(L»„4000)i 


(D„4000)2 


(D„4000)3 


z^ 3.90822 


-11.2489 


10.1242 


-2.60069 


-134.442 


274.564 


-182.224 


39.2976 


z^ 165.269 


-339.815 


227.173 


-48.4853 



Table A. 12. Coefficients for the SDSS i band as a function of redshift (< 0.7) and D„4000 with 

k-corrections derived from the models. 



(Z)n4000)0 


(Z)n4000)i 


(D„4000)2 


(D„4000)3 


24.114 


-53.0651 


38.4687 


-8.97268 


^2 -231.287 


488.825 


-334.771 


74.5591 


z^ 292.203 


-618.56 


423.008 


-93.5608 
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Table A.13. Coefficients for tfie SDSS z band as a function of redshift (< 0.7) and £)„4000 with 

k-corrections derived from the models. 



(L>„4000)0 


(L>„4000)i 


(D„4000)2 


(L»„4000)3 


z^ 14.5525 


-32.5735 


24.1581 


-5.70902 


-203.691 


433.748 


-298.226 


66.4327 


rs' 272.826 


-580.065 


.397. 171 


-87.9759 



Table A. 14. Coefficients for the Johnson-Cousins \J band as a function of redshift (< 0.7) and 
Z)„4000 with k-corrections derived from the models. 



(L>„4000)0 


(Z)n4000)i 


(Z)„4000)2 


(L'„4000)3 


z^ -9.40727 


10.8489 


0.644962 


-1.85733 


z^ -16.6686 


36.2742 


-38.3145 


16.2542 


z^ 113.958 


-246.016 


185.69 


-50.5656 



Table A.15. Coefficients for the Johnson-Cousins B band as a function of redshift (< 0.7) and 
£)„4000 with k-corrections derived from the models. 



(L>„4000)0 


(L»„4000)i 


(D„4000)2 


(D„4000)3 


z^ -30.9335 


56.0725 


-30.7097 


5.99606 


-61.7054 


140.634 


-100.957 


22.3406 


z^ 165.231 


-369.545 


261.328 


-58.2109 



Table A.16. Coefficients for the Johnson-Cousins V band as a function of redshift (< 0.7) and 
-D„4000 with k-corrections derived from the models. 



(Z)n4000)0 


(Z)n4000)i 


(D„4000)2 


(D„4000)3 


z^ 12.8942 


-33.7534 


26.7974 


-6.63462 


z^ -190.187 


396.174 


-262.704 


57.944 


217.218 


-449.526 


298.181 


-65.3664 
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Table A. 17. Coefficients for the Jofinson-Cousins R band as a function of redshift (< 0.7) and 
D„4000 witli k-corrections derived from the models. 



(£»„4000)0 


(L»„4000)i 


(D„4000)2 


(D„4000)3 


12.4246 


-29.1592 


22.4394 


-5.40824 


-175.969 


366.836 


-248.343 


54.7488 


z^ 218.881 


-457.84 


310.701 


-68.0011 



Table A. 18. Coefficients for the Johnson-Cousins / band as a function of redshift (< 0.7) and 

-D„4000 with k-corrcctions derived from the models. 



(D„4000)0 


(Dn4000)i 


(D„4000)2 


(D„4000)3 


z^ 15.6639 


-34.3478 


25.0727 


-5.82947 


-196.21 


412.661 


-280.861 


61.9077 


z^ 256.842 


-540.797 


367.306 


-80.3907 



Table A. 19. Coefficients for the SDSS u band as a function of redshift (< 0.7) and {u — r) with 

k-corrections derived from the models. 





(u - r)" 


(v rY 


(v rf 


[u-rf 


z^ 


2.80197 


-5.1719 


3.01032 


-0.447618 


z2 


-16.0821 


25.7573 


-10.8789 


1.60111 


z^ 


17.9556 


-26.6165 


11.3526 


-1.79009 



Table A.20. Coefficients for the SDSS g band as a function of redshift (< 0.7) and {g — r) with 

k-corrections derived from the models. 





[g-rf 


{g-rY 


{g - rf 


{g - r f 


z^ 


0.184129 


1.00938 


5.40158 


-2.41043 


z^ 


0.112299 


-6.17196 


-7.21778 


4.89703 


z^ 


-2.37339 


15.3331 


-1.95179 


-2.36049 
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Table A. 21. Coefficients for the SDSS r band as a function of redshift (< 0.7) and (g — r) with 

k-corrections derived from the models. 





{g-rf 


{g - r)^ 


{g - r)2 


(g - 




-0.196478 


3.11481 


0.486552 


-1.12639 


^2 


-5.4512 


-6.547 


1.09495 


2.90943 


z^ 


7.81183 


8.41136 


-3.05571 


-2.2226 



Table A. 22. Coefficients for the SDSS i band as a function of redshift (< 0.7) and {g — i) with 

k-corrections derived from the models. 





{g-^f 


{g - iY 


{g - i? 


{g-if 


z^ 


0.0678646 


2.034 


-1.09155 


0.110666 




-7.44388 


4.05582 


0.521655 


-0.136389 


Z3 


10.1914 


-9.5645 


2.07254 


-0.160024 



Table A.23. Coefficients for the SDSS z band as a function of redshift (< 0.7) and (r — z) with 

k-corrections derived from the models. 





ir -.)" 




(r - ^ 


{r zf 


z" 


0.269465 


-0.496672 


1.43921 


-0.220175 




0.935489 


3.91534 


-5.03895 


0.274978 


Z^ 


-2.15467 


-3.71225 


5.19059 


-0.292333 
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3000 4000 5000 6000 ^ 7000 8000 9000 

Wavelength (A) 

Fig. 1.— The 10 CB07 SSP models scaled to unity at 5500 A (top). The ages for the population 
in each model from top to bottom at the blue side of the spectrum are: 5Myr, 25Myr, 100 Myr, 
250 Myr, 500 Myr, 1.0 Gyr, 1.4 Gyr, 2.5 Gyr, 5.0 Gyr, and 10.0 Gyr. The SDSS g and r bands 
and the location of the 4000A-break with respect to these filters for different redshifts (bottom). 
The shaded regions indicate the observed wavelength ranges where we determine L'„4000 for three 
different redshifts. See eq. [2] for the definition of D„4000. 
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Fig. 2. — Integral and differential completeness of SHELS as a function of total i?-band magnitude 
{left axis), and the cumulative number of galaxies with increasing magnitude {right axis). 
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Fig. 3. — Fraction of galaxies where the apertures of the Hectospec fibers {d = I'.'S; solid line) 
and SDSS fibers {d = 3"; dashed line) contain 20% of the total light calculated from the R band 
photometry. 
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'0^?! 0.4' ittf 'O; 200 400 600 

Fig. 4. — Comparison of {g — r) for different sized apertures. ( Top) SDSS 3" fiber {g — r) and 
synthetic Hectospec 1'.'5 fiber {g — r). {Bttom) SDSS Petrosian {g — r) and synthetic Hectospec 
I'.'S fiber (g — r). The solid line traces the median of the differences as a function of redshift. A 
histogram shows the color difference (right). The dotted line indicates the median of the entire 
sample. 
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Fig. 5. — Predictions of the k-corrections for g {left) and r {right) as a function of redshift based 
on the 10 basic CB07 SSP models and two of these models (those with the youngest and oldest 
population) with dust attenuation. The vertical dotted lines indicate the redshift up to which we 
can determine the k-correction from our spectra for these two filters. 
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Fig. 6. — Distribution of D„4000 for SHELS galaxies (top) and the fraction of galaxies with one 
or more emission lines as a function of Z)ri4000 (bottom). Emission-line galaxies have a REW of at 
least 5 A for one or more of the following emission lines: Ha, H/3, [O iii] A4959, [O iii] A5007, or 
[On] AA3726,3728. The thick dashed line indicat es the local ni i nimum of the bimodal distribution 
in D„4000 at L>„4000 = 1.46, very close to that of lWoods et al.l (j20ld ). indicated by vertical dotted 
line at D„4000 = 1.44. The dashed line is a simple tanh fit to the data points to guide the eye. 
The rise at large Dri4000 reflects a contribution by AGN. 
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Fig. 7. — Observer frame {g — r) color as a function of L)„4000 for the 10 CB07 SSP models without 
attenuation (left). The solid lines connect the individual points at equal redshifts. Observer-frame 
{g — r) color as a function of Z)„4000 for the SHELS galaxies (right). The points are color-coded by 
the redshift of the galaxy. The colored solid lines indicate the median of {g — r) binned by redshift 
as a function of -D„4000. Both panels show that {g — r) becomes less sensitive to the age of the 
stellar population of a galaxy toward higher D„4000. 
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Fig. 8. — K-correction derived from the spectra {top and middle) and the models (bottom) as a 
function of redshift where each point is color-coded by the -0^4000 of the galaxy for g (left) and 
r (right). The colored solid lines indicate the median of the k-correction binned by -D„4000 as a 
function of redshift in the top and bottom panels. The black lines in the middle panel compare 
the k-corrections predicted from the model tracks from Fig. [5] with our empirically determined 
k-corrections. 
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Fig. 9. — Assessment of the accuracy of the analytic approximations for the k-corrections in the 
g band as a function of Z)„4000 and redshift. From left to right, the residuals from the surface 
fitting as a function of L'n4000, the distribution of residuals {solid line) and overplotted a Gaussian 
fit to the distribution (dashed line), and the residuals as a function of redshift, with the solid lines 
indicating the median of the difference binned by -0^4000 as a function of redshift. We plot these 
quantities for the k-corrections in the SDSS g band derived from the spectra (top row) and the 
SDSS g band derived from the models (second row). The points are color-coded by the L'„4000 of 
the galaxy. The bottom two rows are the same as the top two rows but for the SDSS r band. 
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Fig. 10. — Assessment of the accuracy of the analytic approximations for the k-corrections in the 
g-band as a function of (g — r) and redshift. From left to right, the residuals from the surface fitting 
as function of {g — r), the distribution of residuals {solid line) and overplotted a Gaussian fit to the 
distribution {dashed line), and the residuals as a function of redshift with the solid lines indicating 
the median of the difference binned by L>„4000 as a function of redshift. We plot these quantities 
for the k-corrections in the SDSS g band derived from the spectra {top row) and the SDSS g band 
derived from the models {second row). The points are color-coded by the Z)„4000 of the galaxy. 
The bottom two rows are the same as the top two rows but for the SDSS r band. 
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Fig. 11. — K-corrections for our galaxies derived from the prescription given bv lCMZlQ based on 
their pegase.2 models for the g (left) and r band (right). Each point is color-coded by the {g — r) 
of the galaxy. 
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Fig. 12. — Comparison of the k-correction derived from our spectra and the prescription given by 



CMZig based on their pegase.2 models for the g (left) and r band (right). Each point is color- 



coded by the {g — r) {top) and the I?n4000 (middle) of the galaxy. The colored solid lines indicate 
the median of the difference binned by Z)„4000 as a function of redshift. The median (solid line) 
and the 68.3% range (dotted lines) around the median for graphs in the top and midd le panel 
(bottom). The vertical dotted line in the left panel indicates the maximum redshift where ICMZIO 
have enough green and blue galaxies to constrain their prescription. 
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Fig. 13. — Test to determine the discrepancy between our k-corrections and those from ICMZIC 
for high Z?„4000 galaxies; difference between the magnitude derived from the spectrum and of 
the model fit for the g (top) and r band (bottom) for the observed (left) and rest frame (right). 
Each point is color-coded by the -D„4000 of the galaxy. The solid lines indicate the median of 
the difference binned by Z?n4000 as a function of redshift. The Figure also shows the influence of 
emission lines on magnitudes; almost all galaxies with low D„4000 and away from the median of 
the difference have emission lines. 
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Fig. 14. — The averaged luminosity-weighted rest-frame spectrum of 84 galaxies with 0.4 < z < 0.6, 
D,i4000 > 1.7 and a difference between the magnitude derived from the spectrum and the fit larger 
than 0.2 {top; see Figure [T3]) . the difference between the average spectrum and fit (middle), and 
the relative difference with respect to the fit (bottom). The dotted line in each panel indicates the 
zero-level. 




Fig. 15. — Test of the accuracy of the model-derived Z?„4000. Relative difference between the 
spectrum- and model-derived Z)„4000 as a function of S/N of the spectrum (top). The solid line 
indicates the median and the dashed line the 68.3 % range around the median. Normalized his- 
tograms of relative difference between the spectrum- and model-derived -D„4000 in different S/N 
ranges {bottom left). The vertical solid and dashed lines indicate the median and 68.3% range 
around the median for each sample. The number in brackets is the number of galaxies in each S /N 
range. Normalized distribution of model-derived Dn^OOO for each of the S/N range {bottom right). 
The percentage in the brackets is the fraction of galaxies with Z)„4000 < 1.46 for each sample. 
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Fig. 16. — Rest frame spectra of five galaxies with different amounts of line emission and (de- 
tectable) emission lines to study the effect of the emission lines on the k-correction. The REW of 
Ho; quoted in each panel is corrected for underlying stellar absorption. Each of these galaxies has 
a redshift of ~ 0.13 — 0.14. The bottom panel shows the transmission curves of the g and r band 
for reference to the rest-frame spectra. 
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Fig. 17. — Influence of the presence of emission lines on the magnitude (top) and k-correction 
{bottom) for the g (left) and r band (right). For each panel the quantity derived from the model is 
subtracted from that derived from the model fit to the data. 
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Fig. 18. — Normalized bandpasses for SDSS ugriz bands {top) and Johnson-Cousins UBVRI bands 
{bottom) as a function of wavelength. The shaded regions indicate the observed wavelength ranges 
where D„4000 would (roughly) be between filters. 
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Fig. 19. — The k-corrections derived from the model fits (top) for the SDSS u band. The colored 
solid lines indicate the median of the k-correction binned by D„4000 as a function of redshift {left) 
and the black lines indicate the k-corrections from models shown in Figure^ (right). Assessment of 
the accuracy of the analytic approximations for the k-corrections (bottom); from left to right: the 
residuals from the surface fitting as a function of Z)„4000, the distribution of residuals (solid line) 
and overplotted a Gaussian fit to the distribution (dashed line), and the residuals as a function of 
redshift where the colored solid lines indicate the median of the residual binned by -0^4000. The 
bottom two rows are the same as the top two rows but for the Johnson-Cousins V band. 



